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A highlighting feature of Majorana bound states in two-dimensional topological superconductors 
is that they gain a phase factor of 7r upon being orbited by a vortex. This work focuses on the vortex 
degree of freedom itself and demonstrates that the change in the Majorana state is accompanied by 
a fermion parity change within the vortex. Such a parity flip is interpreted as a higher dimensional 
analog of the fermion parity pump mechanism in superconducting wires as well as through general 
topological arguments. It is demonstrated in terms of level crossings in three different situations - in 

i) spin-triplet paired superconductors, and in proximity-induced superconducting systems involving 

ii) quantum spin Hall-ferromagnet hybrids and iii) Chern insulators. 


Zero energy Majorana bound states (MBS) are ex¬ 
otic quasiparticles m that support non-local storage 
of quantum information and non-abelian quantum opera¬ 
tions IMS]. They are predicted to appear as topologically 
protected boundary states in several condensed matter 
systems, such as on p- wave superconducting wires [7], 
which can be realized by proximity induced supercon¬ 
ducting strong spin-orbit coupled nanowires in magnetic 
fields TOl - fTB] . ferromagnetic atomic chains on a super¬ 
conductor 03 , quantum spin Hall insulator (QSHI) - su¬ 
perconductor (SC) - ferrormagnet (FM) heterostructure 
junctions [18] and at crystalline defects in two dimen¬ 
sional topological superconductors mm- A highlight¬ 
ing feature is the fermion parity switch induced in such 
Majorana pairs by a phase slip or the encircling of a vor¬ 
tex between them. In this Letter, we explore the “back- 
reaction” of such a switch on the vortex and show that 
a unique feature emerges, namely, the internal states of 
the vortex itself undergo a change in parity. 

More precisely, Majorana fermions are described by 
hermitian operators 7 \ that satisfy the Clifford relation 
{'Yii'fj} = 7i7j + 7 j7» = 25 ij- A pair of MBS encodes a 
two-level system |0) and |1) = c^|0), for c = (71 + *72)/2 
the Dirac fermion operator generated by the MBS pair. 
In a superconducting medium, an electronic quasiparti¬ 
cle acquires a —1 quantum phase when orbiting around 
a quantum vortex of magnetic flux <f> o/2 = hc/ 2 e. As a 
Majorana operator is a linear combination of electronic 
operators, a MBS also picks up a minus sign when a 
well separated flux vortex moves adiabatically around it. 
A paradox now arises from the non-local fractionaliza- 
tion of the electronic degree of freedom into a MBS pair. 
If the flux vortex only encircles one MBS, say 72, the 
Dirac fermion operator is conjugated, c = (71 +*72)/2 
<f = (7i — *72)/2, and the two-level system flips |0) O 
|1). Alternatively, the fermion parity operator, given by 
(—1)^ = *7172, changes sign when 72 —> —72- If the 
pair of MBS and the quantum vortex are well isolated 
from all other low energy modes, their total fermion par¬ 
ity (—l) F +A,ortex cann ot change as tunneling of an elec¬ 
tronic quasiparticle is thermodynamically suppressed by 


the excitation energy gap. The switch of fermion parity 
in the two-level system must therefore be compensated 
by a fermionic excitation at the quantum flux vortex. 
This manifests as a topologically protected level crossing 
among the Caroli-de Gennes-Matricon vortex states EH- 
We refer to this vortex evolution as a fermion parity flip. 



FIG. 1. Numerical results on a 32 x 32 periodic lattice. A 
single level-crossing as a vortex orbits (a) the superconduct¬ 
ing QSH-NI-FM interface eq. §, and (b) the SC trench in 
a Chern insulator eq. (c) and (d) respectively show the 
wavefunction magnitudes of the zero energy MBS pair and 
the vortex state at the crossing for (b). 

It is important to note that in contrast, topological 
phases that support Ising anyons as quantum excita¬ 
tions - such as the Moore-Read fractional quantum Hall 
state [22] and the spinless p + ip superconductor PUHH3] 
- cannot realize vortex parity flips. Each vortex in these 
systems hosts a single MBS and does not associate a lo¬ 
cal fermion parity. Or in more sophisticated language 
the addition of a fermion if to an Ising anyon er does not 


























2 


change its anyon type, i.e. i\> x a = a. Thus, there is no 
measurable change in fermion parity. While an even vor¬ 
tex with even number of Majorana’s carries well-defined 
fermion parity (—lj-Gorte* = i n ^ 1 .. .j 2 n, it does not flip 
when the vortex orbits an Ising anyon as each vortex 
Majorana mode 7, changes sign under the cycle. Vortex 
fermion parity flips are therefore unique in systems where 
MBS bind not to vortices but to static defects. 

In what follows, we first trace the conceptual origin of 
vortex fermion parity flip to the fermion parity pump in 
one-dimensional (ID) p- wave superconductors [71 [Ml[ 25 ] , 
which in turn is the superconducting analog of the Thou- 
less charge pump [ 26 ]. Having established the parity flip 
argument in ID and the associated energy level crossing, 
we explore a range of instances for vortex parity flip in 
two dimensions ( 2 D), each recently proposed as an excit¬ 
ing means of nucleating Majorana bound states. 

As an explicit ID example, the p-superconducting Ki- 
taev wire represented by the lattice Hamiltonian [7 

H — pN = ^ tclcr+i — pc\c r + Aclcl +1 + h.c. ( 1 ) 

r 

is topological and carries zero energy boundary MBS 
when the electron hopping strength |t| is bigger than the 
chemical potential \p \. The low energy states of a su¬ 
perconducting ring with two weak links, one at r = 0 
and the other at r = L/ 2 , are labeled by the two local 
fermion parities (— 1 ) F ° and (—l) Fi / 2 . When the phase of 
the pairing A = |A|e* v winds adiabatically by 27 t along a 
segment, say [0, L/2], there is a level crossing at each of 
the links. This drives the vortex to an excited state after 
a cycle with an extra fermion, which is pumped across 
the bulk although there is a finite bulk pairing gap. 




FIG. 2. Level crossing of a proximity induced superconduct¬ 
ing nanowire when a he/ 2e flux vortex passes across. 

Towards generalizing the fermion parity pump argu¬ 
ment to higher dimension, considering passing a quantum 
flux vortex across a proximity induced superconducting 
nanowire described by Eq. 0. The vortex brings spatial 
variation to the pairing phase A r = |A|e ZV3r , where <p r 
is the polar angle of site r from the vortex core. Fig¬ 
ure [2] shows the level crossing of a 20-site system for 
t = | A| = 2 /r. When the flux vortex crosses the nanowire, 
a Bogoliubov - de Gennes (BdG) state on the nanowire is 
brought down to zero energy with a wavefunction local¬ 
ized at the point where the vortex intersects the wire. At 


the same time the fermion parity of the MBS pair flips. 
This mimics the fermion parity pump because the pairing 
phase winds by 2n within the nanowire segment enclosed 
by the vortex trajectory. After a cycle the bulk nanowire 
is left with a fermionic excitation, which compensates for 
the parity flip of the MBS pair. 
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FIG. 3. (a) Spatial configuration of the d vector around a 
half-quantum vortex (HQV) of a p + ip SC. (b) Zero energy 
Majorana modes of a half-quantum vortex (HQV) and a full 
quantum vortex (FQV). 


Our first proposed 2 D generalization of these argu¬ 
ments is in the context of the spin-triplet chiral super¬ 
conductor having p x + ip y pairing [ 27 H 55 ] celebrated for 
its prospects for supporting non-Abelian vortex Majo¬ 
rana modes [Si- Assuming two spin degrees of freedom 
s =f, I, a half quantum vortex (HQV) [ 36 H 39 ] in such a 
system hosts a zero energy Majorana bound state, 7#, 
at its core. The pairing A ss /cj(cj(, associated with such a 
vortex can be described by A ~ | A\e lv> / 2 {d x — id y , (cr • 
d)a y }e icp / 2 [H], where the Pauli-matrices cr act on spin 
degree of freedom. The phase of the order parameter, <p 
winds by 7r around the vortex and so does d, the direc¬ 
tion of triplet pairing (see figure [ 3 ja)) . The 7r-winding 
of the d vector can be generated by the transformation 
A —► e - l ^/ 4 Ae ze<Tz / 4 , where 9 is the polar angle about 
the HQV. Together with the e* e / 4 phase from the he/ 4 e 
flux, a fermion circling the vortex picks up a minus sign 
for s =f or a trivial phase for s =j.. Thus the zero energy 
Majorana mode 7 h occupies only the spin-up sector. 

Of interest here, the full quantum vortex (FQV), whose 
order parameter winds by 27 t while the d-vector shows 
no winding, contains a spectrum of bound states. Ide¬ 
ally, each of the spin species s =f, J, hosts a zero en¬ 
ergy Majorana yf., in contrast to the HQV, where ef¬ 
fectively the single Majorana state is attributed to one 
of the spin species. Unlike in the HQV, the degener¬ 
ate FQV Majorana modes are fragile towards pertur¬ 
bations and typically hybridize by the effective Hamil¬ 
tonian H = A(— 1 ) F into a level-split ±A energy pair 
with the two members characterized by different par¬ 
ities (— 1 ) F = i'Yp'jp- Here, we specifically consider 
the energy splitting, A, that arises due to a Zeeman 
field, B. In Ref. [ 40 ], which proposed measuring non- 
Abelian statistics via the Aharonov-Casher effect, it was 
shown that A = /iB • d, where /i is the Bohr magne¬ 
ton. For the case where the d vector lies in-plane, in the 
presence of a HQV, the tunnel splitting takes the form 
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H\d\(B x cos(0/2) + B y sm(9/2)). As the FQV encircles 
the HQV as 9 goes from 0 to 27r, the splitting changes in 
magnitude and sign. There is exactly one energy crossing 
when tan(0/2) = —B x /B y and is protected as the FQV 
ground state flips fermion parity. This parity flip can 
be qualitatively understood by seeing both 7 h and 7 ^ 
switch signs after a cycle while jp, which lives on a sepa¬ 
rate spin sector, is unaffected by the process. As a result 
the Dirac operator c = ('y F +i^ F )/2 associated to the low 
energy FQV modes is conjugated to = ( 7 ^ — ij^ F )/2. 

For a similar reason, fermion parity flip occurs in a (p+ 
ip) f x(p — ip) l topological superconductor (TSC). The 
TSC belongs to class Dill as per the Altland-Zirnbauer 
classification m and its topology is protected by time 
reversal symmetry [32] 143] . Upon spin-orbit coupling, 
the two spin-degenerate MBS 7 ^, 7 ^ at a FQV split in 
energy due to the time reversal breaking magnetic flux. 
When the FQV circles around a HQV hosting a single 
MBS in one of the spin sector, there is one (or in general 
an odd number of) level crossing among the vortex states 
in the FQV. 

Fermion parity flips have a general topological origin. 
Their presence is guaranteed by the change of sign of a 
single (or in general odd number of) MBS. In 2D, the 
BdG Hamiltonian U(k, 9,t) that describes the bulk su¬ 
perconducting medium around the MBS varies adiabati- 
cally by the polar angle 9 centered at the MBS as a func¬ 
tion of the adiabatic/temporal parameter t. This class 
D adiabatic cycle has a non-trivial Z 2 topological index 
according to the classification of topological defect [25] . 
The BdG Hamiltonian is topologically equivalent to a 
massive Dirac model 

H{ k, 9, t) = k x T x + k y Ty + mn( 0 , t) + 0 (k 2 ) ( 2 ) 

where T x , Fj,, n(d, t) are mutually anticommuting matri¬ 
ces and r 2 = T 2 = n(0,f ) 2 = 1. The mass parameter 

n lives in the classifying space BO = o(njxO{n) x 
where 0(n) is the orthogonal group and n is related to 
the number of bands in the system 04]. The adiabatic 
evolution defines a map (9,t) —> n(0,f) homotopically 
classified by n 2 (BO) = Z 2 , whose non-zero element char¬ 
acterizes a non-trivial winding and topologically protects 
the fermion parity flip. 

The model © having a slowly varying mass term 
unifies 2D fermion parity flip scenarios in different sys¬ 
tems | |45j . two more of which we now present. The 
first is a proximity induced superconducting (SC) inter¬ 
face between a quantum spin Hall insulator (QSH) [46}- 
149] and a trivial normal insulator (NI). The presence 
of protected zero energy MBS [IE, requires time rever¬ 
sal (TR) breaking and can be facilitated by coating 
an (anti)ferromagnet (FM) along the interface (see fig¬ 


ure 0 b)). We take an 8 -band square lattice model 

H(v,t) = -^QSH-Nl( r ) ® T Z 

+ A x (r,t)r x +Ay(r,t)Ty + h(r)n y (3) 

-ffQSH-Ni(r) = f(sin k x a x + sin k y a v )p x 

+ [m{ r) + e(2 — cos k x — cos k y )\ p z (4) 

where a, and r act on spin, orbital, and Nambu degrees 
of freedom, respectively. The Nambu basis is chosen to 
be (cf^jCj.^, c| M ,-c| iAt ) so that © has a particle-hole 
symmetry 5 = a y T y K, for K the complex conjugation 
operator. Eq.Q describes the QSH-NI interface where 
the mass gap m{ r) changes sign. We assume strong SC 
proximity so that the induced pairing order A = A x +iA y 
is non-vanishing throughout the system. The antiferro- 
magnet couples strongly to a strip neighborhood of the 
interface (see figure 0 where the FM order |/i(r)| out- 
weights the pairing |A| but vanishes elsewhere. 

The QSH-NI interface hosts a gapless ID helical mode 
with opposite spins counter-propagating electrons. The 
helical mode is unstable to TR or charge conservation 
breaking perturbations. Its removal by magnetic field 
was seen in the earliest experiment of QSHI 051 and 
antiferromagnetic (FM) gapped edge was achieved in a 
graphene QSH state m- Moreover, induced SC in QSH 
edge has been observed in HgCdTe quantum wells m- 

FM and SC are competing orders along the interface 
and a FM-SC domain wall - where \h\ — |A| changes sign 
- bounds a protected zero energy MBS. A pair of MBS 
are therefore located at the ends of the ferromagnet in 
figure 0a). The superconducting QSH-NI interface - ex¬ 
cept being TR symmetric and can only be realized holo¬ 
graphically as the edge of a 2D system - can be treated as 
a Kitaev p -wave wire and thus carries protected bound¬ 
ary MBS. When a hc/2e flux vortex passes across the 
superconducting QSH-NI interface, it is akin to traveling 
across a Kitaev p- SC where there is a single protected 
level-crossing among the vortex states. This signals a 
vortex parity flip as the vortex is excited with one extra 
fermion after a cycle. 

We have numerically verified the vortex parity flip phe¬ 
nomenon via its signature level crossing by putting the 
model (0 on a periodic 32 x 32 square lattice (see fig¬ 
ure 0a)). The QSH-NI interface is located along the 
diagonal line and the four sides, which are sandwiched 
between the upper and lower triangular regions with op¬ 
posite insulating mass m. We choose the hopping t = m, 
a uniform pairing strength |A| = 0.5m and the antifer- 
rormagnetic coupling h = 0 . 8 m on a strip over half of 
the QSH-NI interface. To avoid monopole effects 52] . 
we arrange a vortex and an anti-vortex with opposite 
flux, depicted by 0 and 0 in figure 0a). The vortices 
bring spatial and temporal variation to the SC pairing 
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A(r ,t) = |A|e^( r ’*), for 

p Mv,t) = {z - W!(t))(z - w 2 ) , . 

\{z-w 1 (t))(z-w 2 )\ 

where z = x + iy is the complex coordinates for lattice 
point r = (x, y), and wi = Xi + iyi are complex positions 
for the two vortex cores for l = 1,2. The temporal de¬ 
pendence of ([5]) comes from the circular motion of the 
first vortex as it orbits around a MBS when t goes from 
0 to 27r. The second vortex is kept stationary. 

Figure [lja) shows a level-crossing of vortex states and 
confirms the fermion parity flip. At the crossing, a unit 
of fermion is pumped between the vortex and the MBS 
pair. Unlike the p- wave wire case (see figure [ 2 ]) where the 
fermionic excitation is confined along the wire, here the 
excitation stays localized at the vortex as it moves away 
from the QSHI-NI interface until it is brought to the bulk 
bands and hybridizes with the rest of the system. 

As another instance, vortex parity flip can also occur 
on proximity induced superconducting Chern insulators 
(Cl) [53H551 shown in figure [ijb). It can be described by 
the 4-band BdG Hamitlonian on a square lattice 

H =t (sin k x a x + sin k y a y )T z (6) 

+ [m + e(2 — cos k x — cos k y )\a z + A x t x + A y r y 

where a , r again act on spin and Nambu degrees of free¬ 
dom similar to the previous case ©>■ Without the SC 
pairing, ^ describes an insulator with Chern number 
1 when —2e < m < 0. The a z term is a TR breaking 
Zeeman coupling that competes with the induced s-wave 
pairing. We assume the pairing |A| is weaker than the 
insulating mass \m\ so that the 2D system is not in the 
chiral p + ip phase [53] and a full quantum vortex does 
not hold a zero energy MBS. 

The consequence of the bulk Chern invariant is that 
the Cl carries a gapless chiral edge mode that propagates 
in a single direction [531 EH- When two uncoupled Cl’s 
with the same chirality are juxtaposed side by side, the 
interface bounds a pair of counter-propagating electron 
channels with opposite spins ■ This gapless heli¬ 

cal interface can be gapped out by TR breaking backscat- 
tering or U(l) breaking pairing be¬ 

tween the two boundaries. These orders compete and a 
domain wall, where |m| — |A| changes sign, traps a pro¬ 
tected zero energy MBS. This can be realized by insert¬ 
ing a single-layer thick strongly superconducting trench 
in the 2D system, where |A trench | > \m\ (see figure Qb)). 

As with the superconducting QSH-NI interface consid¬ 
ered above, the SC trench in a Cl behaves like the ID 
Kitaev p-wire and flips the fermion parity of a passing 
quantum vortex. We numerically verify this by putting 
the SC-CI model ([b]) on a 32 x 32 periodic lattice. Simi¬ 
lar to the previous case, we arrange a vortex anti-vortex 
pair and consider a circular vortex trajectory around the 
end of the SC trench where a MBS sits. The pairing 


phase A(r,f) = A^, + iA y = |A|e I< ^ r,t ) is also given by 
([5]). Figure [TJb) shows the adiabatic evolution of ener¬ 
gies throughout the cycle. States between ±|m — A| are 
vortex states and are localized at the two vortex cores. 
There is a single level-crossing signifying the parity flip of 
the vortex as it travels across the SC trench. Figure [ljc) 
and (d) shows the localized wavefunctions of the zero en¬ 
ergy MBS pair and the vortex state near the crossing 
respectively. The SC-CI setup is even more prefarable 
than the previous cases in demonstrating the vortex par¬ 
ity flip. Unlike the SC QSH-NI interface, the vortex exci¬ 
tation here stays localized at the vortex core throughout 
the evolution and never hybridizes with the bulk. This 
means that in the absence of accidental fermion poison¬ 
ing, the vortex would carry a different electric charge 
after a complete cycle. 

Finally, we discuss how vortex parity flip has an ana¬ 
log in the recent context of twist defects in topologi¬ 
cal phases with anyonic symmetries [231 I57U5T] . This 
new interpretation is immensely powerful and applies 
even to fractional MBS (or parafermions) [62U65j . Two- 
dimensional s-wave superconductors are fermion parity 
protected quasi-topological phases [53] [37]. They have 
the same topological order [33] as a Z 2 gauge theory [33h 
CD. A quantum vortex of 4 > = hc/2e takes the role of the 
Z 2 flux m, and an excited vortex with an addition BdG 
fermion ip realizes the Z 2 charge e = m xi/j. These quasi¬ 
particles appear in the proximity induced SC QSH-NI in¬ 
terface, and a MBS at the SC-QSH-FM heterostructure 
serves as a dislocation twist defect HHH7] that switches 
e o m when they orbits it. 

SC-CI hybrids on the other hand have a different topo¬ 
logical order. The anyonic content is identical to a U( 1) 2 , 
or equivalently SO( 2)i, theory [61]. A he/ 2e vortex, 
denoted by m, supported by the SC traps a fractional 
charge e* = e/2 (modulo 2e) on the Cl by the Laugh- 
lin argument m and carries semionic statistics. After 
a cycle around a MBS at the end of the SC trench, the 
vortex m is excited with an additional fermion and has 
different charge — e* (mod 2e), and becomes the antipar¬ 
ticle m = m x i/j. A MBS at a SC-CI heterostructure 
can thus be regarded as a twist defect that conjugates 
orbiting quasiparticles. 

In this Letter, we have shown that a highlighting fea¬ 
ture of Majorana defect bound states in two-dimensional 
superconductors, namely the gaining of a phase factor 
of 7T upon being orbited by a vortex, is necessarily ac¬ 
companied by a fermion parity switch in the vortex it¬ 
self. We have presented this scenario in several possible 
geometries relevant to recent theoretical and experimen¬ 
tal explorations of topological systems. In principle, the 
parity flip would be detectable through charge sensitive 
measurements and would constitute not only a signature 
of MBS physics but also a unique parity process in and 
of itself. 
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